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ABSTRACT 

We present spectroscopic observations for a sample of 36 iJer5c/ie[^ -Spire 250-500/im selected 
galaxies (HSGs) at 2 < 2; < 5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Red- 
shifts are confirmed as part of a large redshift survey of iJer^c/ie/- SPlRE-selected sources covering 
~0.93 deg^ in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution 
Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). 
Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts 
at 24 /im or 1.4 GHz using a cross-identification likelihood matching method. Individual source lumi- 
nosities range from log(LiR/L0)=12.5-13.6 (corresponding to star formation rates 500-9000 M© yr~^, 
assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared 
galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and 
infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well 
characterized across large areas of sky in contrast to most z > 2 HyLIRGs in the literature which 
are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, 
we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at 
(7 ± 2) X 10~^ M0 yr~^ h^ Mpc~^ at z ~ 2.5, which is >10% of the estimated total star formation rate 
density (SFRD) of the Universe from optical surveys. The contribution at z ~ 4 has a lower limit of 
3x10"^ M0 yr-^ h^ Mpc"^, > 20% of the estimated total SFRD. This highlights the importance of 
extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, 
even at the earliest epochs. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: infrared — galaxies: star- 
bursts — submillimeter: galaxies 
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1. INTRODUCTION 



Submillimeter galaxies (SMGs, of ten selected by 
850 /im-lmm flux densities ^ 2 mJy; Small, Ivison 



Blain|19971|Hughes et al |1998||Barger et al |1998||Eales 



et al. 1999P are the most intrinsically luminous star- 
burst galaxies that have been identifled to date. SMGs 
are thought to evolve much like local ultraluminous in- 
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frared galaxie s (ULIRGs; [Sanders et al [1988 Sanders & 
Mirabel||1996 | via major mergers. The merger" evolu- 
tionary scenario starts with the cohision of gas-rich disk 
galaxies igniting an intense, short-hved (r ~ 100 Myr) 
phase of gas consumption and dust production via a star- 
burst, fohowed by the formation of a quasar and eventu- 
ahy (1-2 Gyr later) a massive, elliptical galaxy. In con- 
trast to local ULIRGs, SMGs at z ^ 2-3 are much more 
luminous and more massive (in and M//^ ), and some- 
times much larger (Chapman et al. 2004; Biggs & Ivison 
2008J, thus t hey have been dubbed "scaled-up" (Tacconi 



et al. 2008|, provi ding evidence for cosmic downsizmg 
( [Cowie et ar||1996[ ). 

The observation that the most luminous i nfrared 
sources are at the highest redshifts (e.g. GN20; Daddi 



et al. 2009) poses a unique problem for galaxy evolu- 
tion studies. How can these distant ULIRGs be formed 
so quickly after the Big Bang with such high star for- 
mation rates? Their extreme infrared luminosities might 
stem from different evolutionary histories tha n the lo- 
cal ULIRG mergers, i.e. secu lar gas accretion ( Dekel et| 
al. 2009 Dave et al. 2010), but solving the origin of 
infrared- luminous galaxies requires large, uniformly se- 
lected samples of ULIRGs across many epochs. 

Unfortunately, most z > 2 infrared-luminous galaxy 
samples number ~30 galaxies selected in non-uniform, 
biased ways. It is well known that the selection of SMGs 
is severely bias ed, first against ga l axies with warmer dust 
temperatures ([Blain et al. ||2004| [Chapman et al. ||2004[ 



Casey et al. |2doO[[:^(jll[ [Chapman et a l. |:^OlO[ [Magdis et[ 

[2010 ) and second, against galaxies at higher redshifts 



since they are unlikely to have bright radio counterparts 
([Chapman g/. ' 2005) or 24 /im counterparts (Ivison 



et g/. |[2007| [Clements al. |[2008[ ) due to the radio K 



correction and surface brightness dimming. Third, SMGs 
dit z > 2 are often detected in inhomogeneous, serendip- 
itous studies with a range of detection thresholds at dif- 
ferent wavelengths in the far-infrared. Fourth, the spec- 
troscopic follow-up and redshift confirmation of these 
sources is non-uniform; their success rate could relate to 
their FIR properties, e.g. color or single-band flux den- 
sity. The lack of SMGs with confirmed redshifts z ^ 3.5 



several systems at 4 < z < 5.3 (Wang et al. |2007| Wang, 


Barger & Cowie[[2009| |Daddi et al. ||2009| 


Coppin et al. 


1 20091 


Capak et al. 201 1| Smolcic et al. 


[20111 [Walter 


et al. 


2012), however this work has been 


severely lim- 



non-uniform coverage of existing ground-based submil- 
limeter surveys with SCUBA, MAMBO, LABOCA, and 

AzTEC. 

The Herschel Space Observatory ( Pilbratt et al. [2010 ) 
has surveyed ^200 deg^ down to the ^ 5mJy confu- 



sion li mit of Spire (Griffin et al. 



2010) as part of the Herschel Mu' 



2010| [Nguyen et al. 
t i-tiere d Fxtragalac- 



tic^rvey (HerMES; j Oliver et g/~[ [2012D at 250, 350 



and 500 /im. Although high-z infrared galaxies are spa- 
tially rare, SPIRE has mapped much larger sky areas 
than previous submillimeter surveys and thus can detect 
a statistically significant population of z > 2 starbursts 
with a well-characterized selection. 

This paper presents redshifts and spectra for 36 2 < 
z < 5 Herschel SPlRE-selected galaxies identified within 
a large sample of ^ 1600 SPlRE-selected galaxies spec- 



troscopically surveyed over '^Ideg^. In section [2j we 
describe the source selection, biases in 24/im and radio 
samples, and spectroscopic observations. In section [3j we 
present redshift identifications. In section |4j we present 
our results, from derived luminosities, dust tempera- 
tures, the FIR/radio correlation, to composite rest-frame 
ultraviolet and infrared spectra. In section [5j we dis- 
cuss the context of our results by calculating the Her- 
schel SPlRE-selected galaxy (HSG) contribution to the 
cosmic star formation rate density (SFRD), and the im- 
plications for infrared-luminous galaxy evolution in the 
early Universe. In section |6) we conclude. Throughout 



we use a flat ACDM cosmology ( Hinshaw et al. [[2009 ) 
with iJo=71 kms~^ Mpc~^ and r^M=0.27. WTien possi- 
ble, we discuss distance and volume using the general 
unit, Mpc. 

2. SAMPLE AND OBSERVATIONS 
2.1. The Herschel- SPIRE selected galaxy (HSG) sample 
The sources described in this paper were de tected by 



the Spire instrument (Griffin et al. 2010[ ) onboard 



the Herschel Space Observatory as part of the Herschel 



Multi-tiered Extragalactic Survey (HerMES; Oliver et al. 
"2012^. Sources were spectroscopically observed in a large 
redshift survey follow-up prog ram d escribed in detail in a 
parallel paper, [Casey et g/7] ( [2012| ), hereafter C12. The 
results of the redshift survey have been split between 
two papers due to the significant differences in the com- 
prehensive 731 source z < 2 sample, identified through 
[Oil], [OIII], II/3, and Ha emission, than the 36 galaxies 
at z > 2 identified primarily through rest-frame ultravi- 
olet features. We refer the reader to C12 for a detailed 
discussion of our source selection and completeness and 
only briefly summarize those results here. 

Due to the large beamsize of Spire observations (18^^ 
25'^ and 36'^ respectively, at 250, 350, and 500/im), 
counterpart identification and point source pho tometry is 
performed b y extracting flux from Spire maps 



et al. [[2010 ^) at known positions of Spitzer-MUi^'S 24 /im 
and VLA 1.4 GHz sources (see C12 for more details on 
data). This cross-identification prior source extrac tion 
method (called ''X ID" ) is described in detail in ^Rose-[ 



boom et g/.] ( 2010 ) and [Roseboom et al. [ ( |2012[ ). The dis- 
advantage of the XID method is that it relies on Spire- 
bright sources being detectable at 24 /im and/or 1.4 GHz, 
an assumption which is known to sometimes fail at z ^ 3, 
depending on the depth of 24 /im or 1.4 GHz coverage. 
For this reason, our spectroscopic survey was conducted 
in the HerMES coverage areas with the deepest available 
ancillary data in six different legacy fields: Lockman Hole 
North (LHN; a -10h46', S -59°), Cosmic Evolution 
Survey field (COSMOS; S -10hO^ S -2°), Great Ob- 
servatories Origins Deep Survey North field (GOODS-N; 
a -12h36^ 5 -62°), Elais-Nl (ENl; a -16h0^ 5 -54°), 
the UKIDSS Ultra-deep field (UDS; a -2hl9^ S — 5°), 
and the Extended Chandra Deep Field South (CDFS; 
a -3h30^ S —28°). 

Sources were selected for spectroscopic follow-up by 
detection at >3<j significance in at least one of the three 
Spire bands. The absolute flux limit changes field to 
field depending on XID prior source density, but averages 
— 10-12 mJy across the three bands. Throughout the 
rest of this paper we refer to this population as Herschel- 



High-z Herschel-SPlRE Galaxies 



3 



Spire selected galaxies (HSGs) for convenience. Higher 
priority follow-up is given to sources detected in all three 
bands, however the source density of all SPlRE-detected 
sources is low enough such that >98% of all HSGs can be 
surveyed within one spectroscopic mask area, whether it 
be with the Keck I Low Resolution Imaging Spectrome- 
ter (LRIS; covering 5.5' x 7.8') or with the Keck H DEep 
Imaging Multi-Object Spectrograph (DEIMOS; covering 
5'xl6.r). 

2.2. Biases in the HSG Sample 

Although the completeness of the XID source identifi- 
cation teclmiqi^^ < 2 (|Magdis| 
et al |[2QTq1 [Roseboom et al |[2QTq1 [20121 [BeFhermin et| 
[2011 ), the completeness at higher redshifts is difficult 



to estimate since an increasing, non-negligible fraction of 
SPlRE-bright sources drop out at 24 /xm and/or 1.4 GHz 
with increasing redshift. This is also a function of the 
24 jiia and 1.4 GHz depth, which is different field to field. 
Since this paper focuses exclusively on the z > 2 Spire- 
bright population, it is important to emphasize that the 
sample here is incomplete, biased and is not representa- 
tive of all SPlRE-bright galaxies at z > 2. Constraining 
the whole Spire population at z > 2 will require de- 
tailed high-resolution submillimeter follow-up, e.g. from 
ALMA, of a large population of SPlRE-bright systems, 
particularly those that are radio and 24 /im-faint and for 
sources which fail to yield optical spectroscopic identifi- 
cations. 

Note that the purity of the XID technique in coun- 
terpart identification is not guaranteed. In other words, 
XID might be incorrect in its identification of the mul- 
tiwavelen gth counterpart for a Her schel source (as men- 
tioned in Roseboom et al. ]|2Q1Q ). The purity of this 
sample is impossible to gauge without direct far-infrared 
interferometric observations (this does away with confu- 
sion noise, as in Wang et al. ||2Q11 ). However, we do 
note that radio counterpart identitications are more ro- 
bust than 24 jam counterparts due to their source rarity 
and radio's direct scaling with FIR luminosity (Chap-' 
man et al. 2003, |2005|). Half of our sample is radio 
identified, with the other half showing no overall bias or 
trends which would skew our results. 

Another possible bias of the HSG sample is the method 
of spectroscopic targeting. Our Keck observations were 
centered around high-priority 'red' Spire sources (e.g. 
'S'250 < -5350 < 'S'500) which are thought to be the highest- 



redshift Spire sources (e.g. |Cox et al. 1201 1| ). While only 
1-2 'red' targets were chosen per masK, this could skew 
the total redshift distribution higher than if masks were 
laid down arbitrarily on the sky. As C12 describes, and 
as we discuss later in section [5) we measure the impact 
of 'red' sources on the overall redshift distribution to be 
negligible since many of the 'red' targets failed to yield 
spectroscopic identifications. 

2.3. Spectroscopic Observations 

Spectroscopic observations were carried out at the 
W.M. Keck Observatory using the Low Resolution Imag- 
ing Spectrometer (LRIS) on Keck I and the DEep Imag- 
ing Multi-Object Spectrograph (DEIMOS) on Keck II in 
2011 and 2012. LRIS observations were carried out on 
2011-Feb-06, 2012-Jan-26, 2012-Jan-27, and 2012-Feb-27 
with the 400/3400 grism, 560nm dichroic, and primarily 



the 400/8500 grating in the red with central multi-slit 
wavelength 8OOOA for the 5.5^x7.8' mask. This setup 
yields a 1.09 A dispersion in the blue (R--4000) and 0.80A 
dispersion in the red (R'^9000). Integration times var- 
ied from ~2700 to 5600s per mask depending on air- 
mass and weather. DEIMOS observations were carried 
out on 2011-May-28, 2011-May-29, 2011-Nov-28, 2012- 
Feb-16, and 2012-Feb-17. The 600 lines/mm grating and 
72OOA blaze angle was used, resulting in a 0.65A dis- 
persion (R-11000). The GG455 ffiter was used to block 
higher-order light, and typical integration times per mask 
were -2700 to 4800s. The resolution in LRIS red and 
DEIMOS is sufficient to distinguish between a single 
emission line (e.g. Lya) and the [Oil] doublet (rest- 
frame separation of 3A). Data reduction for LRIS was 
done using our own custom-built IDL routines, while we 
used the DEEP2 DEIMOS data reduction pipeline for 
DEIMOS data^^. 

Twenty-five LRIS masks were observed (thirteen un- 
der photometric conditions), and twenty-nine DEIMOS 
masks were observed (sixteen under photometric condi- 
tions), surveying a total of 1594 SPlRE-selected sources 
in 0.93 deg^ (0.43 deg^ observed in photometric condi- 
tions). Of 1594 sources surveyed, 767 have confirmed 
spectroscopic redshifts identified primarily by the [Oil] 
doublet, [OIII], H/3, Ha, [Nil], CaH and CaK absorp- 
tion, the Balmer break, H7, Lya, CIV, GUI], Hell, 
and the Lyman break. The lower redshift sources identi- 
fied through rest-frame optical signatures are discussed 
in C12. Of the 767 confirmed redshifts, 36 are above 
z = 2 and thus comprise the HSG sample discussed in 
this paper. Out of the 36, we categorize 22 as secure 
based on the signal-to- noise of Lya and/or detection of 
multiple spectral features, discussed in the next section. 

Note that the 36 source sample discussed in this pa- 
per was not selected in any special way in Spire color, 
photometric redshift, or optical characteristics. There 
was no special selection imposed which would yield more 
high-redshift identifications, and the selection differs in 
no way from the lower redshift confirmed sources. This 
paper is simply a description of all z > 2 sources con- 
firmed in our large spectroscopic survey. 

The DEIMOS wavelength coverage roughly spans 
45OO-95OOA, whereas LRIS coverage spans 3000-lOOOOA. 
The limited wavelength coverage of DEIMOS results in 
a gap in our redshift coverage between 1.6 < z < 3.2 
which does not occur for LRIS observations. Indeed, 
very few sources observed by DEIMOS are identified 
in that redshift range, with the exception of IHER- 
MESX1.4 J104642.89+585650.0 at z = 2.841 and IHER- 
MESX24J160545.99H-534544.4 at z = 2.555. These are 
both quasars with strong, broad CIV emission. We take 
the DEIMOS redshift desert into account when calcu- 
lating the luminosity function and contribution to the 
cosmic star formation rate density later in Section [5] 
The redshift distribution for all Herschel-SPlRE galaxies 
is shown in Figure [l] The photometric redshift distribu- 
tion shown comes from all HSGs in the COSMOS field, 
not limited to those targeted for spectroscopic observa- 
tions. 

The analysis pipeline used to reduce the DEIMOS data was 
developed at UC Berkeley by Michael Cooper with support from 
NSF grant AST-0071048. 
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SPECTROSCOPIC REDSHIFT 

Fig. 1. — Redshift distribution of HerschelSpiRE selected galax- 
ies from our survey. This paper analyzes the z >2 sample, which 
consists of 36 sources, while the z <2 sample, consisting of 731 
sources, is analyzed in a parallel paper (Casey et al. 2012; C12). 
The distribution of photometric redshifts comes from sources in 
the COSMOS field, where the quality of the photometric redshifts 
is high over a large area; this distribution includes sources not nec- 
essarily in the spectroscopic survey, but follow the same selection 
as spectroscopic targets. The hashed area from 1.6 < z < 3.2 
highlights the DEIMOS redshift range where there is a deficit of 
sources due to spectroscopic incompleteness. 

3. REDSHIFT IDENTIFICATION AND SPECTRA 

Table [l] summarizes the spectroscopic identifications 
and multi-wavelength properties of the 36 SPlRE-selected 
galaxies identified at 2 < z < 5. In each case, the 
signal-to-noise of the identifying feature (Lya emission 
in the majority of cases) is required to be >5 in the two- 
dimensional spectrum. The two-dimensional spectra for 
Lya identified sources is shown in Figure [2] Note that the 
S/N of Lya changes from the two-dimensional spectra to 
one-dimensional extractions of the spectra depending on 
the sources' compactness, the compactness of continuum 
relative to Lya emission, and the observed wavelength of 
the line. For example, one-dimensional extractions with 
Lya longward of ~6000A are prone to contamination by 
the OH forest, thus will be of significantly poorer qual- 
ity in ID than in 2D where OH features are more easily 
distinguished from real lines. 

Initially, we split sources into two categories: sources 
with multiple spectroscopic features (e.g. Lya emis- 
sion, SilV absorption, CIV emission) and single-line 
identifications (where Lya is identified as the only 
emission line). The former group of identifications 
is naturally more secure than the latter. The latter 
source list is dominated by sources at the high red- 
shift end, 3 ^ 2: ^ 5, since multiple features are natu- 
rally more difficult to identify in higher redshift sources. 
However, there are clear cases where the identifica- 
tion of Lya at >5<j is more secure than the iden- 
tification of multiple other features at low signal-to- 
noise (S/N ^ 5cr). An example of a secure single- 
line source is 1 HERMES X24 J160802.63+542638.1 at 
z = 3.415 where Lya is detected at >5<j. An 
example of a less secure multi-feature source is 
1 HERMES X1.4J123622.58+620340.3 at z = 3.579 



where both Lya and CIV are detected at S/N<5cr. In 
this paper, we choose to segregate sources with secure 
spectroscopic identifications from those with tentative 
identifications rather than sources with single- or multi- 
line identifications. 

Sources are categorized as 'tentative' rather than 'se- 
cure' if (a) the signal-to-noise of the Lya feature is 
5 < S/N < 7 as measured in the two-dimensional 
LRIS/DEIMOS spectrum using an appropriately sized, 
adjustable aperture, and if (b) out of five co-authors who 
did a thorough quality ranking of each source's spectrum, 
at least three ranked the source as 'tentative' rather than 
'secure.' Figure [3] show all rest-frame ultraviolet spectra 
for all 36 sources in the sample. 

Tentative sources are marked with a dagger (f 
in Table [l] and denoted as 'tentative' in Figures 
and [3] Out of 36 sources, we classify 22 as secure 
and 14 as tentative. Note that one source, IHER- 
MESX1.4 J100024.00+021210.9 at z = 3.553 is classified 
as tentative not because the Lya emission is of low-S/N, 
but because there is peculiarly strong emission from SilV 
relative to Lya and CIV; this could be an artifact of the 
noise, but since it is particularly unusual, we have cate- 
gorized this source as tentative. 

Photometric redshifts are used, when available, to ver- 
ify spectro scopic re dshifts to within A z / ( 1 + z) ~ 1 / 2 jlT 
bert et al. [2010; Rowan- Robinson et al. 2008; Strazzullo 



et al. ||2010| ,Cardamone et al. ||201Qf 7T:his threshold is 
not strict, out rather used to spot egregious disagree- 
ments between spectroscopy and photometry. These 
photometric redshift catalogs are of varying quality de- 
pending on the depth and number of optical and near- 
infrared bands available, from the limited and shallow 
coverage in ENl to the extensive 30-band coverage in 
COSMOS. The loose constraint of agreement between 
photometric and spectroscopic redshift is based on the 
fact that photometric redshifts are notoriously unreliable 
for dusty starbursts. Even in the deepest multi-band 
fields like COSMOS, ^30% of Spire sources do not have 
photometric redshifts, (e.g. Mobasher et al. [2007 Ilbert 
et al. | [2QTQt. 



Viewmg geometry can have a very strong impact on 
the relative fraction s and wavelengths of escaped ultra- 



violet / o ptical light (Siana et al. 112007, 2010; Scarlata et 



al. [200 9). Heavy obscuration can have draniatic impact 



on photometric redshift reliability in two ways. First, by 
obscuring rest-frame ultraviolet and optical light so much 
that the source has an artificially high photometric red- 
shift (i.e. the source is thought to drop out due to red- 
shifting rather than reddening). Second, dusty sources 
often have differential absorption of Lya photons and 
non-resonantly scattered continuu m photo ns (e.g. the 
"UV chimney" argument, see Neufeld|1991| ), resulting in 
artificially low photometric redshifts. 

Even in the lower redshift samples (C12), sources with 
very confident photometric redshifts from the multi-band 
data of COSMOS and CDFS have photometric redshift 
accuracy of Az/(1 + z)=0.29. This highlights the lack 
of reliability in photometric redshifts of dusty galaxies. 
Therefore, the primary purpose of the photometric red- 
shift restriction in this paper is to verify that the detected 
emission line is in fact Lya and not [Oil] (which would 
correspond to redshifts between and 1). 
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Fig. 2. — Two-dimensional spectra around the identified Lya features for the sources only identified through a single emission line. This 
emission line is not thought to be [Oil] or Mgll in these cases since there is no detection of commonly bright, accompanying emission lines, 
e.g. [Oil], [OIII] and H/3, within the wavelength coverage of LRIS/DEIMOS observations. Note the wide variety of Lya morphologies, 
from very diffuse to very compact, with and without redward continuum. The sources at higher redshifts are observed at higher wavelengths, 
thus likely have one-dimensional extractions contaminated by the OH forest. 
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Figure [2] — continued. 
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Fig. 3. — One dimensional spectra of z >2 Spire sources in observed wavelength shown in the left panels. The flux scaling is arbitrary. The 
wavelengths of rest-frame ultraviolet features are marked by vertical dotted gray lines with noted names. All lines are shown irrespective of 
whether or not they are seen in emission, absorption, and whether or not they are detected. Redshifts are measured off of the Lyo; redshift 
in all cases where Lyo; is detected (a minority of sources is identifled by CIV emission). The Lyo; redshifts are typically redshifted with 
respect to other spectral features, evidence of stellar winds (to be discussed at more length in a future paper). Right panels show best-flt 
far-infrared spectral energy distributions to infrared photometry. The total SED is a solid line while the underlying cold-dust modified 
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Fig. 4. — Spectroscopic redshift against photometric redshift 
for our SPlRE-selected galaxy sample. Sources with identifications 
based on multiple spectroscopic features (emission lines, absorption 
lines, breaks) are shown as black circles, while identifications made 
from a single-feature (i.e. Lyo; emission) are open circles. The 
photometric and spectroscopic redshifts for the full sample of 767 
HSGs (C12) is shown as small gray points. Since the single-line 
source emission lines could have been mis-identified, we also mark 
the corresponding spectroscopic redshifts if [Oil] is assumed and 
find generally poor agreement with the photometric redshift (black 
crosses). The distribution in spectroscopic redshifts for sources 
with no photometric redshifts is shown in gray; there is no obvious 
trend with redshift. 

Fifteen of the 36 sources do not have photometric 
redshifts— even though some of them he in areas of very 
deep ancihary data— a testament to their dusty and op- 
ticahy obscured nature. The sources with fewer identifi- 
able spectral signatures in our spectra are less optically 
luminous than those with multiple features. Those which 
are optically fainter will have less reliable photometric 
redshifts and are more likely to be categorized here as 
tentative. 

3.1. Alternate Redshift Identifications 

Figure [4] plots photometric redshift against spectro- 
scopic redshift where photometric redshifts are available. 
The distribution in spectroscopic redshift of sources with- 
out photometric redshift is the gray histogram. Alternate 
redshifts, assuming the emission lines are actually [Oil] 
instead of Lya, are shown as crosses on Figure [4j The 
[01 1] -implied spectroscopic redshifts are less consistent 
with photometric redshifts than the Lya-implied spec- 
troscopic redshifts. This makes sense in the context of 
the probable enhanced Lya-to-continuum ratios in these 
galaxies. 

While our identified Lya emission lines could also be 
Mgll (2798A), Mgll would be suggestive of a strong 
AGN. With AGN, other features (e.g. GUI] 1909A and 
[Oil]) would also be detectable. Given the expected rel- 
ative strengths of Mgll relative to [Oil], GUI], GIV, 
and Lya, Mgll is unlikely to be identified as the only 
discernible feature. In the case of [Oil], both [OIII] and 
up emission lines and/or GaH & K absorption, or the 
4000A break features are expected. Also at wavelengths 
>6000A, the [Oil] doublet is resolvable with LRIS and 
DEIMOS. In comparison to our lower redshift sample 



(G12), for which there are '^500 [Oil] identified lines 
at the same wavelengths, the sources in this paper have 
none of the above signatures which would point to a mis- 
identification as [Oil]. 

Also note that if the Lya identified lines were actually 
[OIII] or H/3 they would sit at even lower redshifts than 
if [Oil]. Because there are no sufficiently bright emission 
lines between Lya (1216A) and [Oil] (3727A) for star- 
bursts, it is straightforward to segregate between z > 2 
and < z < 1 sources using photometric redshift. 

We inspect the optical images of each source which 
does not have a photometric redshift to judge the plausi- 
bility of our identifications. All three GOODS-N sources 
are i/-band drop-outs. Both IHERMES X1.4 J123622.58 
and IHERMES X1.4J123732.66 are also 6-band drop- 
outs (the two higher redshift sources), and all three have 
optical photometry which is consistent with their spec- 
troscopic identifications. Of the eight LHN sources with- 
out photometric redshifts, IHERMES X1.4 J104557.12 is 
compact and has i ^ 24.5, IHERMES X1.4 J104701.68 
and IHERMES XI. 4 J104709.6 have i ^ 24.3, IHER- 
MES X1.4 J104722.6 and IHERMES X1.4 J104707.7 
have i ^ 25, and IHERMES XI. 4 J104620.4, IHER- 
MES X24 J104642.9, and IHERMES X1.4J104649.9 are 
i-band drop outs. These magnitudes do not rule out 
the possibility that these are low-redshift sources, but 
indicates consistency between our high-z spectroscopic 
identifications and photometry. In GOSMOS, only one 
source, IHERMES X24 J100133.36+023726.9, has no 
photometric redshift since it drops out in all images. 
Elais-Nl sources have much shallower photometric 
coverage than the other fields, hence more sources 
without photometric redshifts. While all sources are 
detectable in wide i-band imaging, multi-band imaging 
is not available across the whole field. None of the 
sources are sufficiently bright or extended at i-band to 
be a convincing z < 1 identiffcation. 

Note that IHERMES X24J161506.65+543846.9, the 
highest redshift source in our sample at z = 4.95 has an 
odd assortment of photometric measurements, dropping 
out in all wavebands (including z-band) except z-band, 
where it has a magnitude of 22.8 (AB). The low photo- 
metric redshift (zp = 1.94 is likely caused by this peculiar 
optical SED, but is also perfectly consistent with Lya in 
z-ban d at z ^ 5, and enhanced Lya-to-continuum ratios 
( Neuf eld|1991 ) . We also note that this source is classiffed 
as tentative. 

4. SAMPLE CHARACTERISTICS 

The importance and context of the 2 < 2: < 5 HSG 
population can only be judged with a basic understand- 
ing on the physical characteristics of the sample. Here we 
measure those physical characteristics, compare them to 
the properties of other galaxy populations, and assess the 
impact of infrared selection biases on our interpretation. 

4.1. SED fits, luminosities and dust temperatures 

The infrared photometry summarized in Table [l] is fft 
with a FIR spectral energy distribution (SED) consisting 
of a coupled single dust temperature modiffed blackbody 
and mid-infrared power law such that 



S{X) = iVbb 



(1 



ghc/>kT _ I 



(1) 
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where S(A) is in units of Jy and T is the galaxy's char- 
acteristic "cold" dust temperature (the dust tempera- 
ture dominating most of the infrared luminosity and dust 
mass). The emissivity index is represented by /3, and Aq 
is the wavelength at which optical depth is unity (here 
fixed at Aq = 200 /im, as described in Conley et aL\ 
2011). The slope of the mid-infrared powerlaw com- 
ponent is given by a, and Ac is the wavelength where 
the gradient of the modified blackbody is equal to a. 
A/'bb and TVpi are the coefficients of the modified black- 
body and power laws respectively. A^pi is a fixed function 
of A/'bb, and T such that the powerlaw and modified 
blackbody are conti nuous at Ac. This SED fitting method 
is described fully in Casey (2012), and is also discussed in 
C12 as applied to the low-redshif t popu l ation. I t is given 
in related forms in Blain et al. (|2002[), Blain, Barnard 
fc Ch apman ( |2003D r |YQunger a7~n[2009p , and [Conley 



,( ,2011^ ). 



To reduce the number of free parameters, we fix a = 
2.0 for sources without Pacs photometry and j3 = 1.5 
for all (the measured values found for local IRAS and 
some distant ULIRGs Casey 2012). These leaves two 
free parameters: A/'bb, which ett'ectively scales with Lir, 
and T, the temperature of the modified blackbody. We 
remeasure dust temperatures for each galaxy by deter- 
mining the wavelength where the flux density peaks and 
convert that to a dust temperature via Wien's Law. This 
provides a more consistent measure of dust temperature 
which can be used in comparisons between SEDs fit with 
alternate techniques, using model templates or direct fits 
jCha ry & Elbaz 2001; Dale & Helou 2002; Blain, Barnard 
^ Ch apman 2003; Siebenmor gejr^FTCr iigel 2007; Draine 
fc Li||2007^ ,Kovacs et al 1,2010 ). We compute infrared 



luminosities by integrating the above best-fit spectral en- 
ergy distribution between 8 and 1000/im. The SED fits 
are shown alongside the sources' optical spectra in Fig- 
ure [3j and their infrared luminosities and dust tempera- 
tures are given in Table [T] 

The luminosities of this sample range from 3.2x10^^- 
6.3xlO^^L0, implying infrared star formation rates 
of 500-9000 M0yr-i (with one outlier at 1.6x10^"^ L©, 
26000 yr~^, whose redshift is tentative). Most galax- 
ies in the sample have SFRs an order of magnitude 
beyond the extreme activity seen in ULIRGs (which 
have SFR ^ 200 -1000 M^t^ yr"^ by the scaling given in 
Kennicutt| |1998|). These starbursts are amongst the 



most extreme star forming galaxies see n in the Universe 



(amongst other HyLIRG populations. Rowan- Robinson 



2000 Bridge et al 



20121 



The conspicuously high star formation rates (e.g. 
above ^lOOOMQyr"^) might lead us to believe that 
AGN contaminate the FIR luminosity or rather, that 
there is potential variation in star formation laws at 
high redshift. There has been some recent discussion 



of whether or not the Kennicutt (1998| scaling between 
IR luminosity and SF R holds under 'ex treme' conditions 
or at high redshifts (Swinbank et al. 12008). Assuming 
a modified IMF would produce more modest SFRs than 
the default Salpeter IMF. While this might change our 
interpretation and change the star formation rates we 
measure here, we use the Kennicutt scaling for SFRs in 
this paper to be consistent with literature work. 

Figure [5] shows the infrared luminosity against dust 
temperature for the HSG sample. There is a noticable 
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Fig. 5. — Infrared luminosity against dust temperature for the 
sample, color coded by redshift interval: 2 < z < 3 (blue), 3 < 
z < 4 (green), and 4 < z < 5 (red). Overplotted are typical lower 
luminosity limit boundaries— a function of dust temperature— at a 
given wavelength, redshift, and flux density limit. A 12mJy flux 
density limit is assumed for the Spire bands (dashed lines) and a 
5 mJy limit is assumed at 850 /im (solid lines). This illustrates how 
half of the SPlRE-selected sample have dust temperatures too warm 
to be 850 /im-detectable. Sources with AGN optical signatures are 
marked with small white dots at their centers. 

absence of very warm sources at lower luminosities. Simi- 
larly, there are very few cold sources at high luminosities. 
This is primarily a consequence of selection effects in the 
SPIRE bands. Warm-dust galaxies are selected against in 
the Spire bands, even at these high-redshifts, due to the 
sensitive variation of infrared flux density measurements 
with dust temperature; the dashed lines illustrate the 
lower luminosity detection limits as a function of dust 
temperature, for a galaxy at z = 2, z = 3, or z = 4. 
This dust-temperature selection bias is even more ex- 
aggerated at 850/imThe luminosity detection limits for 
850/im selection are shown as solid lines on Figure [5j 
nearly bisecting the Spire population so that about half 
would be 850 /im-undetected) . Before Herschel, the cold- 
dust temperature bias of submillimeter observations was 
the focus of many studies looking for the elusive ^Varm- 
dust" SMGs (Blain~a/. ||2QQ4| [Chapman et al ||20Q4 



et 



Casey ani2U09[ |2U1 11 [Chapman et al. ||2U10| [Magdis 
g/. ||2Q1Q[ [Chapin et al. II2Q11I). 

Figure [6] plots HSGs' infrared luminosities against red- 
shift relative to the detection limits of the selection 
wavelengths: the three 250, 350, and 500 /im Spire 
bands, 24 /im and 1.4 GHz. The detection limits in the 
Spire bands depend on the assumed SED shape and 
dust temperature which is not uniform for all Spire- 
selected galaxies. We use the distribution of measured 
dust temperatures in Figure [5] to set limits on the detec- 
tion boundaries in Li^-z on Figure |6] In other words, 
we measure that the mean dust temperature for 250/im- 
selected galaxies to be 40 K, for 350/im-selected galaxies 
as 39 K, and for 500/im-selected galaxies as 37 K. Differ- 
ences on the scale of a few degrees are negligible for these 
illustrative boundary lines in Figure |6j so we adopt a 40 K 
SED for all three selection boundaries. Most of the fields 
which have radio data have a detection limit ^ 40 /iJy, so 
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Fig. 6. — Integrated infrared luminosity (8-1000/im) against 
spectroscopic redshift for our sample (black sources are secure 
and gray sources are tentative, as in Table Luminosities span 
2.4xl0i2_4 0x10^3 ^i^j^ Qj^g outlier at 1.4xlOi^L0. The de- 
tection boundaries of each selection wavelength are also illustrated 
as a function of redshift. The Spire detection boundaries assume 
flux densities >12mJy and an SED dust temperature of 40 K. The^ 
radio boundary assumes that the FIR/radio correlation (Helou, 
ISoifer & Rowan- Robinson 1985 ) holds and a >40 /iJy radio detec- 
tion limit, i'he 24 fim limit assumes a 150 /iJy detection limit along 
with a typical 24/im-LiR scaling (e.g. Le Floc'h et al. 2005). Ra- 
dio sources are marked by a magenta asterisk, while 24 /im sources 
are black. 



we construct the radio detection boundary based on the 



FIR/radio correlation for starburstin^ galaxie s (Helou, 
Soifer fc Rowan- Robinsonj] 19851 Condon 1992) with giR 
evolving as in Ivison et al. ( |2010p . The 24 /im detection 
boundary is the least certain as it scales with Ljr; many 
recent works note up to ~1 dex disagreement between 
extrapolated 24 /xm infrared luminosities and direct mea- 
surem ents (e.g. |Le Floc'h et al ||2QQ5l [2QQ9| jElbaz et al 



2Q11[) . This tells us that, although we can use the 24 //m 



detection boundary shown in Figure |6] as a rough guide, 
it should not be concerning that 24 /im-selected sources 
fall below the line by ^ 0.3 dex. 

Although very few of these sources have existing 
850 /im data, this work (especially Figure [5| suggests 
that half of the SPlRE-selected population would be un- 
detectable in the original ScuBA 850 /im surveys even at 
these redshifts. Specifically, we estimate that only 31/36 
sources would have Ss5o > 2 mJy (86%), while only 21/36 
sources have 6'85o > 5mJy (58%). The statistics of the 
temperature-bias selection effect are discussed in more 
detail in C12, as they relate to the lower redshift popu- 
lation where the statistics are more robust. 

4.2. FIR/radio correlation 

For the 23 adio-detected sources, we investigate the 
FIR/radio correlation for starburst galaxies ( ,Helou„ 
Soifer fc Rowan-Robinson 1985 Condon 1992). Mea- 
suring the FIR/radio correlation in this sample is use- 
ful for checking that our sample is roughly consistent in 
FIR/radio given expectation from our measured redshifts 
(i.e. it is another reassurance on tentative identifications 
in particular). This correlation is measured via the ratio 



qiR log 



1.02 X IO^^Lfir 
-log(l X 10-^^S'i.4(l 



cm 



(2) 



where LpiR is the far-infrared luminosity measured in 
the range 40-120 /im given in L©, Di^ is the luminosity 
distance in cm. Si. 4 is the 1.4 GHz fiux density in /iJy, 
and a is the synchrotron slope, here set to 0.75 (Ibar 



et al. 



2010 



Ivison et al. 2010) and defined such that 



Twenty-three of the 36 galaxies in our sample are radio 
detected (64%), and their measured giR ranges from 0.7 
to 2.1 with mean value (<3'ir)= 1.58 ±0.35. Note that the 
infrared luminosity component in qi^i is I/fir(40 — 120), 
not Lir(8 - 1000). If Lir(8 - 1000) is used instead, lu- 
minosities and q increase by 0.40±0.15dex and the scat- 
ter in q grows; the increased scatter is caused by the 
contribution from the mid-infrared fiux to Ljr. Only 
one galaxy in our sample is 'radio-loud' and indicative of 
AGN (the source at 2; = 3.579); this is consistent with the 
observed CIV emission in its rest-UV spectrum. Since 
we observe the rest of the radio-detected sample to agree 
within uncertainti es with previous nieasures of qm(z) in 



previous samples ([Ivis on et al. |2010[|Kovacs et al. [2010 

Magnelli et al. 2012) lends additional credence to our 
redshift identifications. While each of the literature sam- 
ples have different measures for q (ranging 1.3-2.2), all 
with uncertainties on the order of ~0. 15-0.25, the overall 
trend of an evolving (71R is consistent between samples. 

4.3. Composite Ultraviolet Spectra 

Since the signal-to-noise on individual galaxy spectra 
shown in Figure [s] is quite low for most sources (except 
in the detection of Lya), we construct a composite rest- 
frame ultraviolet spectrum which can serve two purposes: 
it validates bulk redshift identification by way of detect- 
ing lower S/N spectral features around Lya, and it begins 
to shed light on the intrinsic rest-frame ultraviolet emis- 
sion properties of extremely infrared luminous starbursts. 
Unfortunately, larger samples are necessary to perform 
the latter analysis; in this work, our primary goal is to 
help validate our redshift identifications through cross- 
correlation to the composite spectrum for sources with 
only single-line identifications (e.g. Lya). 

It is clear from Figure [3] that these infrared-starbursts 
exhibit a wide range in spectral properties, from Ly- 
man Break Galaxy (LBG) spectra, quasar spectra, star- 
burst spectra, to those with very steep to very shallow 
UV slopes. We classify sources as LBGs if they exhibit 
a steep cutoff in continuum fiux at rest-frame 12 16 A, 
as quasars if they exhibit broad, high-ionization emis- 
sion lines, and starbursts as narrow-emission line galax- 
ies whose emission line luminosity is more significant 
than continuum luminosity. Many sources are dust ob- 
scured and thus are noisy except for the detection of 
Lya emission. The construction of a composite spec- 
trum serves as a sanity check on the redshifts. The 
detection of lower S/N features in a composite spec- 
trum does not directly confirm that every galaxy added 
in is correctly identified, but it does indicate that most 
of them are. Since co-adding spectra from higher S/N 
sources would wash out the low S/N features of the 
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other galaxies (in a sample of 36), we exclude the fol- 
lowing sources from the composites (excluded on the 
basis of detection of non-Lya spectral features at >3<j 
significance, including continuum and CIV): IHERMES 
X24 J033136.96 -275510.9, IHERMES X24 J095917.28 
+021300.4, IHERMES X24 J095948.00 +024140.7, 
IHERMES X1.4 J100008.64 +022043.1, IHERMES 
X24 J100020.16 +021725.2, IHERMES X24 J100036.00 
+021127.6, IHERMES X1.4 J100024.00 +021210.9, 
IHERMES X1.4 J100111.52 +022841.3, IHERMES 
X24 J100146.56 +024035.6, IHERMES X1.4 J104620.40 
+585933.4, IHERMES X1.4 J104636.00 +585650.0, 
IHERMES X24 J160545.99 +534544.4, and IHERMES 
X1.4 J033151.94 -275326.9. In other words, the com- 
posite spectra are only made up of Lya single-line detec- 
tions, those with low ~2-3(T CIV detections, and those 
without. 

Two different composites are constructed based on the 
detection or non-detection of CIV emission at this low 
S/N level. Since most starbursts are expected to show 
absorption in CIV, the co-addition of sources with and 
without CIV might easily yield a null result and no ab- 
sorption or emission. Six sources are co-added in the 
CIV emission composite (these are the remaining sources 
for which CIV emission is detected, as indicated in Ta- 
ble [l] attributed to an AGN). The remaining 17 sources 
are co-added to form the composite without CIV emis- 
sion. Each composite is constructed by scaling the flux 
of each galaxy to an arbitrary fixed mean value in the 
wavelength range 1330-1400A. This wavelength range 
is chosen for its proximity to Lya and absence of spec- 
tral signatures. The two composite rest-frame ultraviolet 
spectra are shown in Figure [T] 

Both composites show a very high-S/N Lya line, Lya 
break, SilV absorption, and either CIV emission o r ab- 
sorp tion. The composite LBG spectrum from S hapley et 
\al. I ([2003) is shown for comparison, although the rela- 
tive sample sizes should be contrasted. The CIV emis- 
sion composite, by design, consists of galaxies with non- 
negligible AGN emission; as a result, the width of the 
Lya line is broader in the CIV emission composite than 
in the CIV absorption composite. 

We test for consistency between the composite 
spectra and individual source spectra through cross- 
correlation in the off-Lya wavelength regions (note 
that for each galaxy's spectrum, we measure the 
cross-correlation with a composite excluding that 
source). This provides an indication of sources which 
might be contaminating the composite rather than 
boosting its signal-to-noise. All individual sources 
in the CIV composite have correlations > 0.4 at 
AA = offsets. The individual source spectra making 
up the other composite are of lower signal-to-noise, 
and have cross-correlations ranging 0.2-0.7. However, 
four sources (1 HERMES X24 J095830.24+015633.2, 
1 HERMES X1.4 J104707.69+585149.1, 
1 HERMES X1.4 J104649.92+590039.6 and 
1 HERMES X1.4 J123536.28+623019.9) have almost 
no correlation with the composite (0.05-0.15) which 
is caused by no continuum detection (since only the 
off-Lya spectra are considered in the cross-correlation 
test). Three of these sources are considered 'tentative' 
in their spectroscopic identifications in Table [l] The 
fourth source, 1 HERMES X24 J095830.24+015633.2, is 



secure as judged by the quality of the Lya detection and 
the inconsistency of this line being incorrectly identified 
as [Oil]. 

4.4. Spectral Signatures of AGN 

Since many of the sources in our sample have clear 
AGN features in their optical spectra, one might think 
that the infrared luminosities are contaminated by signif- 
icant AGN heating rather than starburst heating. Typ- 
ically, the presence of an AGN warms dust to tempera- 
tures ^ 100-200 K. In this sample the majority of galax- 
ies have dust temperatures ^ 70 K. Furthermore, two 
QSOs and the 13 sources with CIV detections have dust 
temperatures in the 30-50 K range, perfectly consistent 
with star formation dominated infrared emission. While 
there is still potential for AGN contribution to Lir, the 
lack of correlation with dust temperature indicates that 
the effect is small ( ^ 25%, the nominal contribution 
of mid-infrared powerlaw emission to Lir, Casey|[2QT2 ). 
This is consistent with prior measures o f SMGs with 
AG N ranging ^15- 25% (IS winbank et al |[2Q04, Pope et 
al 1120081 [Menendez -Delme stre et aL ||2009| |Laird et aL 
20101 ICoppin et aF|2010|) . Note that one study, Alexan- 
der et aL. (|2005|), could b 
this work ( 



DC interpreted to disagree with 
ncLmg --75% of SMGs have AGN), however 
a minority of the sources in that data have AGN which 
dominate the sources' bolometric luminosity. 

We can draw some basic conclusions from the CIV 
emission and spectral types in our sample to infer the 
overall AGN content of z > 2 HSGs. For sources of suf- 
ficient signal-to-noise (>5cr in continuum), we can assess 
AGN spectral signatures source by source. Of the 20 
galaxies which meet this S/N cut, there are three Lyman 
Break Galaxies (LBGs), three quasars, seven starbursts 
with AGN (e.g. CIV emission), and seven "pure" star- 
bursts (see Table [l] for details). Of the remaining single- 
line identifications, none show AGN signatures. Of the 
36 sources, 10 have AGN signatures, three of which are 
obvious quasars. Although very qualitative, this anal- 
ysis implies an AGN fraction of ~ 25% for the z > 2 
HSG sample. From the composite spectra (which was 
constructed from only the lower luminosity sources) , our 
statistics agree by construction; in other words, 6 out 
of 23 sources had CIV emission, or ~26%. In a series 
of detailed studies on the multiwavelength propertie s of 
850 / im-sel ected SMGs , A lexander et al. (2005), Pope etj 



oTl (|2008|) and iMenendez-Delmestre et oTTpOOPI also 

— r ^ _i w-^. 



measure AGN fractions ~ 25% for similarly luminous 
z ^ 2 starbursts. 

Interestingly, the sources exhibiting AGN signatures in 
the optical do not show hotter dust temperature SEDs in 
the infrared. One might expect higher dust temperatures 
in the infrared with the presence of an AGN heating the 
surrounding material to temperatures ~100-200K, ex- 
ceeding normal heating from star formation, ~30-50K. 
The sources with AGN signatures are marked with small 
white dots in Figure [5] The observation that the AGN 
does not seem to have significant impact on infrared lu- 
minosity or dust temperature is not surprising if you con- 
sider that the star formation activity is at least an order 
of magnitude more luminous. 

Although the selections of the SMG population and 
the HSG population differ, finding 1/4 with AGN might 
suggest that HSGs are similar in most ways to SMGs 
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Fig. 7. — A comparison of the rest-frame ultraviolet co mposite spectra for He rschel-selected galaxies (HSGs) with the composite of 
Lyman Break Galaxies (LBGs; top spectrum) compiled in [Shapley et q/.~| (1 2003') from 811 individual galaxies. The composite spectra 
of HSGs are split into those with CIV emission (middle spectrum) and (JiV absorption (bottom spectrum) and exclude all sources with 
high-S/N (>3cr) spectral features which are not Lya emission. There are six galaxies in the CIV emission stack and 17 galaxies in the CIV 
absorption stack, both having a median redshift of {z) =3.6. Lyo; and the Lyman break, SilV absorption and CIV emission/absorption 
are detected in both composites. More observations of similar sources is needed to enhance the S/N of the composites, which exhibit a 
wide range of spectral properties (as seen here and in Figure [sj. The spectra are arbitrarily offset in flux with marked zero-points. 



without any enhanced AGN activity, despite shghtly 
warmer overah dust temperatures and brighter mid- 
infrared fluxes in comparison (described in the next sec- 
tion). However, further detailed work on these samples 
is needed before any conclusion is drawn as to the evolu- 
tionary nature of these iJer5c/ie/-selected galaxies relative 
to classic 850 /im-selected SMGs. 

4.5. Composite Infrared Spectra 

Figure Is] combines all Spire, Spitzer dind Pacs (where 
available jinfrared photometry for all galaxies in our sam- 
ple from rest-frame ^40 to 150/im and radio data. This 
includes observations at 24/im, 70/im, 100/im, 160/im, 
250/im, 350/im, and 500/im. At these redshifts, the 
Spire bands probe the Wien-side of the thermal dust 
emission peak. Mean SEDs are fit using the modified 
blackbody plus power law method described by Equa- 
tion [T] for photometric data which are scaled to the mean 
infrared luminosity of the sample, (i^m) = 1.8 x 10^^ L©, 
and then separately, scaled to the mean radio flux den- 
sity of the sample, (5'i.4) = 95 /iJy (or rest-frame Si. 4 = 
283 /iJy assuming a = 0.75). While all 36 galaxies are 
used in the former SED fit (top panel of Figure [8| , the 
radio-scaled SED fit only has contributions from radio- 
detected galaxies. 

There are two notable aspects of these mean SED fits 
seen in Figure [Sj the first is the difference between the 
observed 24 /im flux densities relative to predictions from 
a 850 /im-selected SMG template spectrum ( jPope et al 
|2008), and second is the difference in dust temperatures 
Detween luminosity- and radio-scaled SEDs. The issue 
of the discrepancy of mid-infrared emission relates to 
the ongoing discussion of suppressi on of PAH em ission 
in infrared starburst galaxies (e.g. Elbaz et al. I2Q11| ), 
where it is suggested that most normal galaxies nave a 
fixed ratio (= IR8) between 7.7/im emission (or ^8/im 



emission) and total integrated infrared luminosity, Ljr 
and that sources with enhanced infrare d emission are 
c alled infrared starbursts. However, as [Hainline et aL\ 



[ ( 2009[ ) point out, the mid-infrared portion of the spec- 
trum does not lend itself to simple interpretation in terms 
of what is or is not AGN dominated or starburst domi- 
nated. SMGs, known to be extreme starbursts, do in fact 
have enha nced infrared emission relative to PAH strength 



e.g. see ,Pope et al. ||2008| [Menendez-Delmestre et al. 



2009), however in this sample— which is on average more 



distant than most 850 /im-selected SMGs— we see mid- 
infrared flux densities ~2-5 times the SMG expectation, 
more consistent with the measured IR^ value for most 
"main sequence" galaxies. What does this suggest about 
the Spire sample's evolutionary histories? Is PAH emis- 
sion simply not suppressed in these distant starbursts, 
or could the 24 /im "excess" be due to AGN heating? 
Or could these high redshift infrared luminous galaxies 
be "main sequence" secularly evolving galaxies? While 
this might be a selection bias based on 24/im or 1.4 GHz 
detectability, the radio-selected galaxies (bottom panel) 
still show a mid-infrared excess above expectation from 
the SMG template. 

The differences in dust temperatures between the two 
mean SED fits (42±3K vs. 54±4K for luminosity-scaled 
and radio-scaled, respectively) is traceable to a radio se- 
lection bias. For two galaxies with similar Spire flux 
densities, one with a warm temperature (~60 K) and one 
with a cold temperature (~30K), the galaxy with the 
warm temperature is going to have a much higher inte- 
grated infrared luminosity and therefore much brighter 
1.4 GHz detection at these redshifts. Therefore, when we 
consider just radio-detected galaxies, the average dust 
temperature increases due to the exclusion of cold, non- 
detectable galaxies. 

As discussed earlier, we estimate that 16-43% of these 
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Fig. 8. — The mean infrared and radio SEDs for the sample re- 
normahzed to the mean IR luminosity of the sample, l.SxlO-*^^ Lq 
(top) and to mean radio flux density, 95 fiJy (bottom). The mean 
normalized flux density in log(A) =0.1 bins is shown as black 
squares, from rest-frame ^4-150/im. Best-flt SEDs (solid black 
lines)are generated as described in the text according to Equa- 
tion m with flxed /5 =1.5. They comprise a cold-dust modifled 
blackbody flt (dashed line) and a mid-infrared power law represen- 
tative of warm dust emission. Radio synchrotron emission is added 
onto this best-flt infrared SED by assuming the FIR/radio corre- 
lation holds with a synchrotron slope of a =0.75. The composite 
SMG SED described in Pope et al. (2008) is shown as a dashed 
line. Both the luminosity-scaled SED and radio-scaled SED ap- 
pear to have a 24 /xm excess relative to the SMG expectation. The 
radio-scaled SED has a hotter characteristic dust temperature than 
the luminosity-scaled SED (both of which are uncertain by ~3K), 
likely driven by the bias against colder-dust galaxies of similar flux 
densities (less likely to be radio detected). 

galaxies would be formally undetected at 850 /im at < 2- 
5 mJy. In other words, 16-43% of HSGs at 2: > 2 are 
consistent with the Submillimeter Faint Radio Galaxy 
(SFRG, formerly optically fai nt, ^^OF RG") selectio n and 
not S MG selection (Chapman et al. |[2QQ4^ ^Casey et al. 



20091 IChapman et al. ||2U1U| [Magdis a/. ||2U1U| [Casey 
et oL ||2011|?| ). The composite infrared SEDs from Fig- 



ure [5] support this conclusion, since the range of observed 
850 /im flux densities from the best-flt SEDs is in the 
1 < 5*850 < 10 mJy range, not as luminous at long wave- 
lengths as the 850 /im-selected composite (Pope et al. 
[2008, ). 

5. SFRD IMPLICATIONS AND DISCUSSION 

To place these z > 2 SPlRE-selected galaxies in context 
with other high— z infrared galaxies and lower luminos- 
ity galaxies, we estimate their contri bution to the cosmic 
star formation rate density (SFRD; pMadau et al. ||1996| 
Madau, P ozzetti & Dickinson 1998; Hopkins ^ BeacornJ 



the build-up of stellar mass in the Universe over a range 
of epochs. At lower redshifts z ^ 1, ultraluminous in- 
frared gal axies are very rare a nd contribute little to the 
SFRD (Sanders et al || 2003 [) but towards z ^ 1, the 
importance oTULIRGs grows, and it is estimated that 
LIRGs and ULIRGs (Ljr >10^^ Lq) could contribute as 
m uch as ^1/2 of the to tal SFRD (see results from Spitzer; 

3 Floc'h et al l|2005| and work in C12). At z > 1, the 
contribution of infrared-luminous sources is much more 
diflicult to measure, lim ited by small numbe rs of SMGs 
( Chapman et al. [2005 Wardlow et al. |20"TT ) or complex 
selection biases or extrapolations from themid-infrared 



(ICaputi et al. ||2007[ iMagnehi et al. f 
20111). The 2 < z < 5 galaxies in th: 



2011 Capak et al. 



2006). The SFRD contribution allows a direct compar- 



ison of the importance of infrared-luminous galaxies to 



IS paper provide a 
unique sample to make this measurement, due to their 
well characterized selection over a relatively large sky 
area, ~ 1 deg^. 

To arrive at SFRD estimates, we flrst compute the in- 
frared luminosity function using a 1/Fmax method, where 
each source is associated with the maximum volume in 
which it could be detected at its given luminosity, Lir 
( [Schmi dt 1968). The number density of sources with lu- 
minosity between L and L + AL is given as ^z{L)AL = 
l/(Vi{L) X Ci) in units of Mpc~^ logL"^. Here q is 
a completeness estimator which corrects for sample in- 
completeness at the selection wavelength, e.g. in this 
case, at 250-500/im. C12 presents a detailed discus- 
sion of this completeness factor as a function of selection 
wavelength flux density which varies fleld to fleld (based 
on the prior source catalog depths). Sources with flux 
densities >15mJy will be more than 90% complete (e.g. 
Ci > 0.9) in ah flelds. 

The completeness estimator's effect on the integrated 
SFRD is small compared to the uncertainty in the lumi- 
nosity function itself from small number counts. Note, 
however, that this correction only pertains to one of the 
many sources of incompleteness of this sample; it is far 
more difficult to quantify and correct for incompleteness 
with respect to sources missing from the prior catalog 
list (e.g. 1.4 GHz or 24 /im faint), which is particularly 
a problem at z > 2, as well as spectroscopic incomplete- 
ness, i.e. the galaxies which have no detectable emission 
lines or are too optically obscured to be identified. 

Vi{L) is the maximum volume in which source i could 
reside and still be detectable by our survey. Since the de- 
tection limits of Spire alter between field catalogs, the 
detection limit determining the maximum volume is cal- 
culated source by source. This luminosity detection limit 
is determined much like the curves in Figure |6] For ex- 
ample, a source might have its highest S/N at 350/im, 
then its luminosity detection limit, thus maximal red- 
shift limit ^max, is determined by setting a 3<j detection 
threshold at 350/im where sigma is the local confusion 
plus instrumental noise in the 350/im map. This Zmax 
limit is then found across the entire survey area probed 
to determine accessible volume. The assumed dust tem- 
perature is that measured for the given source (since dust 
temperature does impact the steepness of the luminosity 
limit with redshift). 

We split the luminosity function into two redshift bins: 
2.0 < z < 3.2 and 3.2 < z < 5.0 with 15 sources 
in the former and 20 in the latter. Since LRIS and 
DEIMOS have different wavelength coverage, DEIMOS 
observations suffer from a redshift desert from 1.6 < 
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Fig. 9. — The estimated luminosity function for 2< z <5 
f/er5c/ie/-SPlRE selected galaxies compared to the luminosity func- 
tion of 850 /im-selected SMGs at z ?^2.5 (Chapman et ||2005l 
[Wardlow et al. 2011). Up arrows denote the fact that this survey 
IS spectroscopzcally incomplete and that the incompleteness is not 
well quantified at z > 2. Numbers next to each point indicate how 
many galaxies from our sample contribute to that luminosity bin; 
the numbers are comparable to those in the z ~2 SMG samples. 
The luminosity function for the whole sample at 3.2< z <5.0 is 
shown in salmon, while the high confidence (h.c.) identifications' 
luminosity function is shown in dark red. 

z < 3.2 that LRIS observations do not, so we split the 
sample at z = 3.2 and only compute the density us- 
ing LRIS observations between z = 2.0 and z = 3.2. 
This excludes two sources from the calculation (IHER- 
MESX1.4 J104636.00+585650.0 at z =2.841 and IHER- 
MESX24 J160545.99+534544.4 at z =2.555) which were 
both surveyed with DEIMOS and detected in the red- 
shift desert due to strong CIV emission caused by the 
presence of a quasar. Since high-redshift z ^ 2 sources 
are only detectable on the masks observed in the best 
weather conditions, the effective area probed by LRIS for 
this calculation is 0.13 deg^ over the range 2.0 < z < 5.0 
and for DEIMOS 0.30 deg^ over the range 3.2 < z < 
5.0. The resulting spectroscopically incomplete luminos- 
ity functions are shown in Figure [9] 

Since some of our identifications are less confident than 
others (e.g. those marked with a f in Table [T]) we 
also compute the luminosity function excluding tenta- 
tive identifications. The result is seen in Figure [9] while 
the 2.0 < z < 3.2 luminosity function remains the same 
(differing only by two sources), the 3.2 < z < 5.0 lumi- 
nosity function drops by 11 sources (salmon vs dark red 
LFs). 

The SFRD is then the luminosity weighted integral of 
the luminosity function, or the raw summation of the 
luminosity (converted to SFR) over accessible volume: 
(7.O±2.O)xlO-^M0yr-ih3Mpc-3 at 2.0 < z < 3.2 and 
(5.5±0.6) X 10-^ Mo yr-i h^ Mpc-^ at 3.2 < z < 5.0 (full 
sample) and (2^^) x 10"^ M© yr"^ h^ Mpc"^ at 3.2 < 
z < 5.0 (high confidence sample), shown in Figure [lO] 
against other comparison populations. These points are 
lower limits since they do not include any sources which 
might be excluded from the prior catalogs at 24 /im or 
1.4 GHz, which is speculated to be a non-negligible frac- 
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Fig. 10. — Star formation rate density of HerschelSPlRE s elected 
gal axies (blac k points) relative to 850 /xm-selected SMGs ( Cliap-| 

fn et al. '2005 Ward low et al. 2011) and 1.2 mm-selected MMUs 
)seboom et al. 2U12 1). I'he co impilation of SFRD measurements 
n (Hopkins &; Beacom 2006) is shown as a gray band, which is 
largely drawn from optical or rest-frame ultraviolet-selected galaxy 
populations corrected for dust extinction. The luminosity limits of 
integration are io(i2-4-i3.4) Lq at z ~ 2.5 and io(i2-8-i3.6) 
z ~ 4; the SFRD from sources sitting in the luminosity where the 
two redshift bins overlap (10(12.8—13.2) j^^^ -g g^own in green. The 
SFRD measurements for Herschel-selected galaxies at z <2 are 
shown as gray points. The sharp drop in the SFRD of Herschel- 
selected samples at z ~2 is caused by the redshifting of the SED 
peak such that more infrared-luminous galaxies are Spire "drop- 
outs" and that only the warmest, most-luminous >10-'^'^ Lq sys- 
tems are detectable with Spire at 2; > 2. Note that we observe an 
increase in the infrared-luminous contribution to the SFRD from 
z ~ 2.6 to z ~ 4. 
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tion ( > 20%) at z > 2 {e.g 

Small etal. in prep). Note also that the luminosity lim- 
its of the two redshift bins differ: the z ^ 2.5 bin covers 
^q12.4-13.2Lq while the z - 4 bin covers lO^^-^-^^-^ L©. 
To assess luminosity evolution from z ^ 2 — 5, we com- 
pute the SFRD contributions in the overlapping lumi- 
nosity range of both redshift bins, IQ^^-^^^^-^ L©, shown 
as green points in Figure [lOj 

Although these measurements of the Herschel contri- 
bution to the SFRD at 2 < z < 5 are lower limits due to 
our survey's incompleteness, the effects of gravitational 
lensing and clustering could lead to an overestimation. 
Are these effects significant in this sample? For the for- 
mer we use the conditional lensing probability as a func- 
tion of Spire flux den sity (Francesco De Bernardis, pri- 
vate communication, [Wardlow et al 2012[ ). For Spire 



flux densities £'500 < 80mJy the distribution in number 
counts is dominated by a Schecter function rather than 
the flat-sloped source counts at > 80 mJy; galaxies with 
SbQQ > 80 mJy have a high probability of being lensed 
by factors > 2 while this model predicts a mean lensing 
factor for this sample of (/i) < 1.05, which changes neg- 
ligibly between z = 2 and z = 5. 

The brightest source in our sample with S'250 = 
73.1 mJy, IHERMES X24 J161506.65+543846.9, has the 
highest probability of being lensed (its expected lensing 
factor is (/i) = 1.2). It is the highest redshift source in 
our sample at z = 4.952. Due to its extreme luminos- 
ity compared to the rest of the sample and its tentative 
spectroscopic identification, we exclude this source from 
the star formation rate density (SFRD) measurement. 

To assess the impact of clustering on the SFRD mea- 
surement, we need a good grasp of the spatial density 
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of z > 2 sources on our slit-masks and the possibility 
of biased placement of slit-masks around high-z clusters. 
The former can be gauged by the number of high— 2: con- 
firmations per slit-mask; the 36 sources of this sample 
are distributed across 12 LRIS masks and 14 DEIMOS 
masks, with an additional 5 LRIS masks and 3 DEIMOS 
masks without any high-z sources; this averages to 0-2 
galaxies per mask without any mask having more than 2 
sources. Since the masks were distributed randomly with 
respect to one another in each field and none of the masks 
were close together, this demonstrates that these sources 
are indeed randomly distributed over the surveyed area, 
0.93 deg^ for the whole survey. 

The possibility exists that there is an intrinsic bias of 
the placement of our slit-masks such that more high-z 
sources are observed than elsewhere. As explained in 
detail in C12, masks were placed around high-priority 
targets which were 'red' in their Spire colors (e.g. 
5*250 < ^350 < ^500 all with S ^ 15 mJy) and thought 
to be high redshift sources. Of the 36 confirmed z > 2 
sources, seven (~ 19%) were originally high-priority tar- 
gets. However, an additional 44 high-priority targets 
were identified at z < 2 and 78 were unidentified. Rela- 
tive to the number of high-priority sources targeted, we 
measure 5±5% as identified at 2: > 2, 34±5% at z < 2 
and 60±4% unidentified. The same statistics for lower 
priority targets are 2±1% at z > 2, 48±1% at z < 2, 
and 50±1% unidentified. Within uncertainties, the pro- 
portion of sources identified at z > 2 are the same be- 
tween low-priority and high-priority targets, indicating 
no bias or advantage in targeting 'red' sources more than 
any other significant Herschel-SPlRE source. This implies 
that no clustering correction on the measured SFRD is 
necessary. 

The lower limits to the SFRD set by Spire sources tells 
us that the early Universe potentially had a very sub- 
stantial amount of star formation in short-lived, intense 
>lOOOM0yr~^ bursts as opposed to slow-progressing 
moderate levels of star formation. The contribution from 
Spire to the SFRD at these epochs is at least compara- 
ble to the contribution measured from longer wavelength- 
selected galaxies, like the 850 /im-sele cted SMGs (Chap- 
man et al 120051 [Wardlow et al |201 1|) or 1.2 mm-selected 
MMCs ( IkosiBoom et al. ||2012p . Th is is made more in- 
teresting D3rT[ie""olDservi!I^^ the populations (SMG 
and HSG) only overlap by 21 out of 36 galaxies (58%). 
Further work aimed at confirming redshifts oi z >2 Spire 
sources, particularly those without radio or 24 jiia coun- 
terparts, is needed to constrain these lower limits into 
real measurements so the importance of ^ 10^^ Lq ac- 
tivity in the first few Gyr of the Universe is understood. 

6. CONCLUSIONS 

The identification of submillimeter galaxies at early 
epochs in the Universe's history is the key to understand- 
ing the limits of star formation and galaxy evolution on 
short timescales. This paper has presented new obser- 
vations of 36 HerschelSviKE selected starburst galaxies 
between 2 < z < 5, taken from a large Keck spectro- 
scopic survey of 1594 SPlRE-selected galaxies covering 
0.93 deg^ 

We present the following conclusions: 

• Our sample of 36 iJer5c/ie/-selected galaxies (HSGs) 



constitute some of the brightest, most extreme in- 
frared starburst galaxies in the Universe. Span- 
ning 2 < z < 5, our sample has a mean luminosity 
(Lir) = 1.8 X 10^3 Lo (SFR^31OOM0yr-i). 

• These 2 < z < 5IISGs have a well character- 
ized selection across 6 legacy fields and 0.93 deg^; 
galaxies must be >3<j significant in one of the three 
Herschel-SPlRE bands and also be detected in deep 
24 /im and/or 1.4 GHz survey coverage. Although 
it misses 24/im or 1.4 GHz high-z dropouts, the se- 
lection is identical to low- 2: HSG selection. Sources 
at z > 2 comprise 5% of all galaxies selected via this 
method; although inefficient for finding high-z in- 
frared galaxies, the selection is easily reproducible 
and well suited for volume density estimates. 

• Our sample show a wide range of rest-frame ultra- 
violet spectral features: some galaxies classifiable 
as quasars, some as LBGs, and most as starbursts 
with a wide range of dust extinctions/reddening. 
The heterogeneous nature of their spectra pro- 
vide additional evidence that the infrared-luminous 
stage might exist during a period when the host 
galaxy is rapidly evolving. 

• The radio-detected subset of our sample (23/36) 
follow the FIR/radio correlation consistent with 
moderate evolution in qi^i from previous work ( |Ivi-] 
son et a/.l2010D . 



• We construct composite rest-frame ultraviolet 
spectra and rest-frame infrared spectral energy dis- 
tributions to assess some aggregate properties of 
HSGs. In the rest-frame UV, we determine that 
25% of HSGs exhibit CIVemission (a signature of 
AGN). In the infrared, HSGs exhibit a 24/im ex- 
cess relative to SMGs of similar Ljr; without mid- 
IR spectra, it is impossible to know whether this 
is due to enhanced PAH emission (similar to "nor- 
mal" galaxies) or AGN emission. The dust temper- 
atures of radio-selected samples are warmer than 
those of the full sample. 

• Our spectroscopic survey is incomplete due to se- 
lection bias at 24/im and 1.4 GHz, as well as spec- 
troscopic incompleteness caused by heavy dust ob- 
scuration in the rest-frame UV. Therefore, we are 
able to place lower limits on the contribution of 
2 < 2; < 5 HSGs to the cosmic star formation 
rate density, which is >7xl0~^ Mq yr~^ h^ Mpc~^ 
dit z ^ 2.6 and >3x 10"^ yr-^ h^ Mpc"^ at 
z ^ 4, corresponding to >10% and >20% of the 
best-estimates of the total SFRD at their respec- 
tive epochs. 

This work highlights the importance of extremely lumi- 
nous FIR-bright galaxies to the build-up of stellar mass, 
particularly at early times in the Universe's history. Fur- 
ther work on constraining completeness and the parent 
population of infrared-luminous galaxies at 2; > 2 is 
needed to understand the role that short-lived starbursts 
have in the context of galaxy evolution and formation. 
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